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Abstract 



ff^ . Assuming that the space-time is close to isotropic in the sense that the shear 

^ I parameter is small and that the maximal velocity of the particles is bounded, we 

have been able to show that for non-diagonal Bianchi I-symmetric spacetimes with 
collisionless matter the asymptotic behaviour at late times is close to the special case 
of dust. We also have been able to show that all the Kasner exponents converge to 

r^ ' 3 and an asymptotic expression for the induced metric has been obtained. The key 

^^ ■ was a bootstrap argument. 

The sign conventions of [18] are used. In particular, we use metric signature — h + 
+ and geonietrized units, i.e. the gravitational constant G and the speed of light c are 
set equal to one. Also the Einstein summation convention that repeated indices are to be 
summed over is used. Latin indices run from one to three. C will be an arbitrary constant 
and e will denote a small and strictly positive constant. They both may appear several 
times in different equations or inequalities without being the same constant. A dot above 
a letter will denote a derivative with respect to the cosmological (Gaussian) time t. 

1 Introduction 

For both mathematical and physical reasons one would like to be able to understand 
anisotropic and inhomogeneous cosmological models in general. As a first step spatially 
homogeneous spacetimes can be studied. For the case where the matter model is a perfect 



fluid a lot of results have been already obtained which are summarized in [22], see also [B] 
for a critical discussion. Usually in observational cosmology it is assumed that there exists 
an Era which is "matter-dominated" where the matter model is a perfect fluid with zero 
pressure, i.e. the dust model. A kinetic description via collisionless matter enables to study 
the stability of this model in the following sense. Suppose we have an expanding universe 
where the particles have certain velocity dispersion. One might think that due to the ex- 
pansion the velocity dispersion will decay. That this is true has been shown in [12], [IS] 
and [TTj for locally rotationally symmetric (LRS) models in the cases of Bianchi I, II and 
III. These results have been generalized (in a different direction) to LRS Bianchi IX and 
in a context which also goes beyond collisionless matter in [2] and [1]. 
The LRS models can be diagonalized in a suitable frame where they then stay diagonal 
if one makes some extra assumptions on the distribution function. A natural question 
is, what happens if the model is not diagonal? In this paper we will treat the late time 
dynamics of the non-diagonal Bianchi I case assuming small data. That this makes sense 
was established in |13| where it was shown that geodesic completeness holds for the general 
Bianchi I-Vlasov case. In a sense which will be specified later, we assume that the universe 
is close to isotropic and that the velocity dispersion of the particles is bounded. Then we 
conclude that the universe will isotropize and have a dust-like behaviour asymptotically. In 
fact these two properties are intimately linked in the proof, something which is not expected 
to happen in other Bianchi types, except in the case of a positive cosmological constant 
where it has been shown [9] that isotropization and asymptotic dust-like behaviour occurs 
in all Bianchi models except Bianchi IX and without the LRS assumption. 
We hope to be able to extend this result in the near future a) to other Bianchi types, where 
the spacetime does not (necessarily) isotropize, but some other solution still may act as an 
'attractor' b) to remove the small data assumption(s). 



2 Bianchi I spacetimes with colhsionless matter 

A Bianchi spacetime is defined to be a spatially homogeneous spacetime whose isometry 
group possesses a three-dimensional subgroup G that acts simply transitively on the space- 
like orbits. They can be classified by the structure constants of the Lie algebra associated 
to the Lie group. We will only consider the simplest case where the structure constants 
vanish, i.e. the case of Bianchi I with the abelian group of translations in R^ as the Lie 
group, where the metric has the following form using Gauss coordinates: 

^g = -dt^ + gabit)dx''dx\ (1) 

We will use the 3+1 decomposition of the Einstein equations as made in [18]. We use 
Gauss coordinates, which implies that the lapse function is the identity and the shift vector 
vanishes, so comparing our metric with (2.28) of [IB] we have that a = 1 and (3'^ = 0. The 



only non-trivial Christoffel symbols are the following (See (2.44)- (2.49) of 

(2) 
(3) 
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where p, ja and T^^ are the energy density, matter current and energy momentum tensor 
respectively. In the 3+1 formulation the second fundamental form kab is used and rewriting 
its definition as in (2.29) of [l^ we have: 

Qab = -2kab. (4) 

Using the Einstein equations as in (2.34) of [IB] and the fact that gab does not depend on 
the spatial variables, we have: 

kab = Hkab - 2kackl - M{Sab " -^^Qab tl S) - ^TlpQab (5) 

where we have used the notations tr 5" = g°'^Sab, H = ti k. It has been assumed that the 
cosmological constant vanishes. From the constraint equations (2.26)-(2.27) of |18] : 

- kabk"'' + H'^ = IQixp (6) 

and we see that no matter current is present since the "spatial" Christoffel symbols vanish: 

Toa = 0. (7) 

From ([6]) we see that H never vanishes except for the Minkowski spacetime. We exclude 
this case and assume now without loss of generality that H < for all time. If this is not 
true for a given solution it may be arranged by doing the transformation t i— )■ — t. 

For the matter model we will take the point of view of kinetic theory. This means 
that we have a collection of particles (in a cosmological context the particles are galaxies 
or clusters of galaxies) which are described statistically by a non-negative distribution 
function /(x°,p") which is the density of particles at a given spacetime point with given 
four-momentum. We will assume that all the particles have equal mass (one can relax this 
condition if necessary, see |2]). We want that our matter model is compatible with our 
symmetry assumption, so we will also assume that / does not depend on x". In addition 
to that we will assume that there are no collisions between the particles. In this case the 
distribution function satisfies the Vlasov equation (See (3.38) of 



f + 2*J/^ = 0. (8) 



where f is defined on the set determined by the equation 

called the mass shell. The energy momentum tensor is (compare with (3.37) of |18]): 

P= fit,p){m^ + gabp'^p^t^idet g)-^dp (9) 

Sab = f f{t,p)paPb{m^ + gabP^'pThdetgy^dp (10) 

Toa= [ f{t,p)pa{detg)-^dp (11) 



Here p := {p^,p^,p^) and dp := dp^dp^dp^. For this kind of matter all the energy conditions 
hold. In particular p > ti S > 0. Our system of equations consists of the equations (|1])- 
(ITT]) . For a given Bianchi I geometry the Vlasov equation can be solved explicitly with the 
result that if / is expressed in terms of the covariant components pi then it is independent 
of time. This has the consequence that if to is some fixed time and fo{pi) = f{t,Pi) then 
we can express the non-trivial components of the energy momentum tensor as follows: 

P= fo{pi){m^ + g'''^PcPd)Hdetg)~^dpidp2dp3 



2.1 Other equations and new variables 

A useful relation concerns the determinant of the metric ((2.30) of |18]): 

^[log(det^)] = -2i7 (12) 



Taking the trace of the mixed version of the second fundamental form (2.36 of |18]): 

H = H^ + AntrS-Unp (13) 

With ([6]) one can eliminate the energy density and flT3l) reads: 

if =i(i/2 + 3A;^,r')+47rtr5 (14) 

We can decompose the second fundamental form introducing Cab as the trace-free part: 

H , , 

Kb = CFab + —gab (15) 

Then 

kabk''' = aabcr'"' + ^ (16) 
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and ( I14p takes the following form: 



H = ]^H^ + Ka,a'^'' + AntiS (17) 

It is obvious that 

H > -H^ (18) 

From the constraint equation ([6]) and (IT7|) it also follows that 

H KH"^ (19) 

From (fT8|) and ( TT9|) we conclude that 

-(t + C)~^ >if(t) > -2(t + C)-^ (20) 

Let us use the following notation: 

F='-^. (21) 

This quantity is related to the so called shear parameter, which is bounded by the cosmic 
microwave background radiation and is a dimensionless measure of the anisotropy of the 
Universe (See chapter 5.2.2 of |22])- Using this notation and with the help of (j3]), (^ and 
(IT7|) we have 

F = H[F{l-^F-^^trS)- ^-^S^t^^'] (22) 



From the constraint equation 



- - levr^ 
3 H^ 



^ = o - 16vr^ (23) 



In general we see that < F < |. 

2.2 Special cases 

2.2.1 Vacuum 

We see that in the vacuum case F = | and F = 0. Actually in this case one can write 
explicitly the solution known as Kasner solution: 



where the constants pi are called the Kasner exponents which satisfy the two Kasner 
relations given by: 

Pl+V2+Vi = 1 

iPif + {P2r + iPsf = 1 
The mean curvature in this case is 

H = -r^ 

Let Xi be the eigenvalues of kij with respect to Qij, i.e., the solutions of 

det(A;} - XS}) = (24) 

We define 

P^ = ^ (25) 

as the generalized Kasner exponents. They satisfy the first but not in general not the 
second Kasner relation. 

2.2.2 Flat Friedmann with dust 

A special case of Bianchi I is the Friedmann model with fiat spatial geometry. The metric 
in this case is: 

^g = -de + a\dx^ + dy"" + dz^) 

where a is the scale factor. In this case aab = which means that F = and also that 
F = 0. 

Moreover in the dust case Sab = (the Einstein-dust system can be thought of as singular 
case of the Einstein- Vlasov system, see chapter 3.4 of |18| for more information) and we 
can solve f lT7|) obtaining: 

H = -2t-^ 

We can also write down the explicit solution of the metric in this case: 

4^ = -dt"^ + t^dx^ + dy"^ + dz^). 

This is also called the Einstein-de Sitter model. 



2.2.3 Bianchi I: the dust case (with small data) 

Let us look at the dust case not necessarily isotropic. The general solution is known 
and one can see from the solution (11-1.12) of [12] that the spacetime will isotropize. 
Nevertheless we will analyze this case with care, since we will show that the general case 
behaves asymptotically like it assuming small data. In the dust case: 

F = HF{1 - -F) < (26) 

We are interested in the asymptotic behaviour at late times and will assume that F is 
small, i.e.: F < e. Then it follows from (IT] 



Integration leads to 



H={\ + ^-F)H'<i^- + e)H' (27) 



H < 



t + C 
Using this inequality in fl26l) and integrating: 

F = 0{t-'^+'). (25 

We can put the equality (l27j) in the following form: 



where Id is 



H-' = -H{to)-' + ^(t - to) + Id (29) 



r lF{s)ds (30) 



Id 

J to 

Now we will use the fact that we can choose freely the time origin setting to = — 2f/(to)^^. 
Note that for the general Bianchi I symmetric Einstein- Vlasov-system we know that H 
takes all values in the range (— oo, 0) (Lemma 2.1 of jT3]) We then obtain 

H + 2t-^ < AlDt-'^ (31) 

Using ( 128|) it is clear that Id < C'(to). Note that with our time origin choice C = in (120|) . 



so 

H > -2r^ (32) 

Our result for H is the following: 

H = -2r^[l + 0{r^)] (33) 



This equation in (12B|) leads after integration to: 

F = 0{t-^) (34) 

From (132|) and ( 134|) we can conclude that: 

aaba^' = 0{t-^) (35) 

Now from f lT5]) we see that the eigenvalues fl2^ of the second fundamental form with 
respect to the induced metric are also the solutions of 

det(ai - [A - ^-H]S'^ 

Let us define the eigenvalues of ctjj with respect to gij by Aj, we have that: 

A,- = A; H 

3 

Note that Sj(Aj)^ = aab'^"'^- From (132|) and ( !35|) we can see that the spacetime isotropizes 
at late times, in the sense that 

where pi are the generalized Kasner exponents. 
Now using (IT2|) and (132|) we have: 

det^ = 0(t^) 

4 

So the hope would be to show that \gab t~3\ can be bounded by a constant. Let us define 

_ _4 

gab = t -igab 



^ab ^ ^+|^afe_ 



We have then with (HI) that: 



gab = --{2t ^ + H)gab-2t ta^^ 
Making similar computations as in ^ we arrive at: 

Il^a6(t)|| < \\-gab{to)\\+ [ [l\2s~' + H{s)\+2{a,ba''\s))-^]\\gab{s)\\ds 
and with Gronwall's inequality we obtain: 

Wgabim < \\-gab{to)\\^w{ I [\\2s-^ + H{s)\ + 2{aaba^\s))^ds} < C (36) 



Therefore gab is bounded for all t > Iq. The same holds for g"-^ by similar computations. 
Thus: 

\t~^^gab\ < C 

\t+lg''^\ < C 
From this we can conclude that: 

Ikafell < cri 



or 



O-ab = 0{t 



Looking again at the derivative of gab and putting the facts which have been obtained 
together, we see that: 

tab = 0{t-^) 

This is enough to conclude that: 

9ab = t+l[gab + 0{t-^)] 

where Qab is the limit of gab as t goes to infinity. 

3 Bianchi I: the general case 

For the general case the basic equations are (IT7|) . a modified version of (|5]) where 7 is a 
small and positive quantity which is introduced for technical reasons and 



d.{-H-') = l + lF + 4.'j^ (37) 

i-(rV'') = rVl^(2t"' + H)- 7t-i] + 2t-^+ta'^^ (38) 

F = H[F{1 -\F-^^ivS)- ^^a.a'^l (39) 

We have a number (different from zero) of particles at possibly different momenta and 
we will define P as the supremum of the absolute value of these momenta at a given time 
t: 

P{t) = suplbl = (rtaP6)^|/(t,p) ^ 0} (40) 



Theorem 1. Consider any C°° solution of the Einstein- Vlasov system with Bianchi I- 
symmetry and with C°° initial data. Assume that F(to) and Pito) are sufficiently small. 
Then at late times one can make the following estimates: 

H{t) = -2ri(i + 0(r^)) (41) 

P(t) = 0{t~l+') (42) 

F{t) = 0{t~^) (43) 

Remark (H2]) implies that asymptoticaUy there is a dust-like behaviour (see (H7|) ). 
Proof We will use a bootstrap argument (see chapter 10.3 of [12] for more information). 
Let us look at the interval [tQ,ti). Our bootstrap assumptions are the following: 

F{t)<A{l+ty^ (44) 

P{t) < B{1 + t)-T2 . (45) 

where A and B are positive constants which we can choose as small as we want. 
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(l37|) we obtain: 
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(46) 



's)]ds 

where the inequality comes from (j33j) and ([6]). Consider now an orthonormal frame and 
denote the components of the spatial part of the energy-momentum tensor in this frame 
by Sah- The components can be bounded by 



so we have that 



from which follows that: 



Sab < P\t)p (47) 



^ < 3P2 (48) 



Ig<^^{A + B^) 



Now returning back to (I46p and setting as in the dust case tg = ^2i/(to) ^ we have besides 

H > -2r^ (49) 
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that 

H + 2t-^ < 18{A + B^)t-^ (50) 

2. Estimate of P 

Now we will use the second equation f l38|) to obtain an estimate for the metric. One 
can show that in the sense of quadratic forms the following is true: 

Then with the definitions of g"-^ and F and the estimates flH|) and f H^ we have: : 

2^-i+t^ab ^ ^^-^+1 ^^^^^ab^i ^ab ^ 2F^H\t-^g''^ < AA^{1 + tyh-H-^g"^ 
Now from the last inequality of the first step (l50l) equation (138|) leads to: 

^(r^r") < r^r1l2(A + i?2)to + 4A^(1 + to)-3 - 7]ri = -r^t-^r't-' 
A and i? have to be chosen in such a way that r] is positive, then 

|(t"^r') < (51) 

from which follows that 

r' = oif) 

or 



^'*'' = 0{r 



Now from ( 15T]) we have: 



7+1 



t-^+-^g'^\t)<t,^^-'g'^\to) 



Using the fact that pa is constant along the geodesies we can conclude that: 

Pit) < to-^+^P{to)t^~^^ < Bt^-^ (52) 

since we can choose P(to) and B independently as small as we want. In order to improve 
(H5|) 7 has to be smaller then |. Using the notation C = i ^^^ ^^^^ inequality can be 
expressed as 

P{t) < Bt~l+< (53) 

11 



where C < n- 

3. Estimate of F 

Until now we have an estimate for H and for P in the interval [to,ti). Now we have 
to improve the estimate for F coming from the bootstrap assumption. The desired esti- 
mate is F{ti) < A(l + ti)~^"'"^. If this is the case (case 1) the bootstrap argument will work 
and there is nothing more to do. 




Bootstrap assumption 
Desired estimate 

1 

Case I 

Let us suppose now the opposite, that F{ti) > A{1 + ti)"^"*"^. Then define ^2 as the 
smallest number not smaller than to with the property that F{ti) > A{1 + ti)^'^'^'^. In this 
case, we have to distinguish between the case that ^2 = ^o (Case Ila) and ^2 > ^o (Case 
lib). Let us look now at ( 15^ . in particular at the terms in square brackets. By using 



3„ 3,, ,3 3,, ,3 

-F<-A(l + t)-2<-A(l + to)-^<- 



where 5 is a positive small constant. Using 



and 



H^ - 2 - 2 ^ 

With the Cauchy-Schwarz inequality, ([6]) 
interval [t2, ti]: 



167r ^ < F2- — 



ty-^ <-B\i + to 



7 S 

-6 < _ 

- 4 



fH7|) and supposing that F > A{1 + 1)""^^^ in the 

5 



H^ 



< V3FF~^p^ < VsA-^B^n + ty^~^F < -F 

P 4 



Note that although A 2 may be a big quantity, since A and B are independent we can 
make B smaller to "correct" this. Using flSU]) and the last three inequalities in (15^ leads 
to: 

F < {-2 + £){!- 6) Ft"^ 
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where e = 18{A + B'^)tQ ^. Now setting ^ = e + 25 — e6 we end up with: 

which means that 

F{h) < F{t,)t^-kr'+^ 
Case Ila. In this case ^2 = ^Oj so ( I54p means 

F{h) < F{to)to'-k,-'^^ < At,-'^^ 
since we can choose F(to) as small as we want. So it follows that in this case: 

Ati~^+' < F(ti) < Ati~^+^. 
Situation which is schematically depicted in the following figure. 




to=t2 



Bootstrap assumption 

Obtained estimate 
Desired estimate 
1 



Case Ila 



(54) 



Case lib. In the case lib we can use the fact that by continuity -^(^2) < A{1 + ^2) ^"""^ 
holds and then 

Fih) < Ail + t2)-'+^t2'-«tr'+« < A{i + hy-ki-^^^ 

The e here is also a quantity which we can choose as small as we want and then it follows 
that in this case 

F(ti)<A(i + tor«tr'+« 

We can choose e to be smaller then ^, so 

F{ti) < Ah-^+^ 
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Bootstrap assumption 

Obtained estimate 
Desired estimate 
1 



Case lib 



4- Results of the Bootstrap argument 



We have arrived at the statement that at least for a small interval [to, ^i) and assuming 
we obtain the estimates: 



Fit) < At-^+^ 
P{t) < Bt-^i+^ 



(55) 
(56) 



Thus both assumptions (144|) and (145|) have been improved and using a bootstrap argument 
we know that the estimates obtained are valid for ti = oo assuming small data. (150|) can 
be expressed as: 



H = -2t~^{l + 0{t~^)) 
and is then also valid for the whole interval. 



(57) 



5. Improving the estimate of F 

We want to improve (155|) , but before that we need an inequality in the other direction. 
From (139|) we have: 

F> H{F + V3F^P^) 

Implementing now (IS^ -( 15B|) in this inequality, in particular having in mind the 'B'-dependence 
of the estimate of P (l52l): 



F>HF + HC{to)BH 



2^-l+e 



Using now (H9]) : 



F> -2t-^F-C{to)BH 



9 10 ^t 
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from which follows 



fF{t) > tlF{to) - f C{to)Bh-l+'ds 

Jtn 



'to 

and 

F{t) > t-\F{to)tl - B'C{to)to-^+^ + C{to)BH-^+^) > t-\F{to)tl - B^ C {t,)t,-^+^) 
Fito) is strictly positive. Choosing now B small enough we have the following estimate: 

Fit) > C{to)r^ (58 

Putting the last term of fl39|) with the help of (l56l) and (158|) in the following manner: 



|167r^^| < CF{F--^P^) = FO(ri+^) 



we can improve (155|) implementing (l55l) - fl57j) in (139|) with the result: 

F = 0{r^) (59) 

With these results we obtain the following theorem. 

Theorem 2. Consider the same assumptions as in the previous theorem. Then 

Pi = l + 0{t~') (60) 

and 

9ab = t^HGab + 0{t~')] (61) 

^"'' = t-t[^"^ + 0(t-2)] (62) 

where Qab and Q"-^ are independent oft. 

Proof The conclusions made in the dust case after ( 134|) only depend on the estimate of F 
and H and apply directly to the general case. 



4 Conclusions and Outlook 

The result concerning the asymptotics of the metric implies "asymptotic freezing" in the 
expanding direction in the following sense. Consider the metric (4.12) of [4], then by 
comparison with (16T]) one sees that the off-diagonal degrees of freedom ni,n2 and n^ tend 
to constants and thus are not important for the dynamics. See [^ for the importance 
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of asymptotic freezing in the "other" direction which has been studied, namely the initial 
singularity, and [^ for consequences of that in a quantum version. 

As already mentioned in the introduction there exist several results concerning the Bianchi 
I-symmetric Einstein- Vlasov system (without cosmological constant). Almost all of them 
assume additional symmetries namely the LRS and the reflection symmetry condition 
(see jT^ for a precise definition of these symmetries). However there exist some results 
where only the reflection symmetry is assumed. Concerning the expanding direction there 
is theorem 5.4 of [13]. Our theorems can be seen as a generalization of that theorem since 
we obtain the same the result, but a) we also obtain how fast the expressions converge b) 
we obtain an asymptotic expression for the spatial metric c) we do not assume any of the 
additional symmetries mentioned. However we used a different kind of restriction namely 
the small data assumptions. 

In any case we think it interesting to study the non-diagonal case because although this 
time there was not an essential difference in the result with respect to the diagonal case, 
we do not expect that this will always be the case, especially when analyzing the initial 
singularity. In [5] the possible dynamical behaviour towards the past has been determined 
assuming only the reflection symmetry and already there surprising new features like the 
existence of heteroclinic networks arose. 

The situation of several (tilted) fluids leads naturally to the non-diagonal case. The Bianchi 
I-symmetry implies the absence of a matter current such that a single tilted fluid is not 
compatible with this assumption. However there can be several tilted fluids such that the 
total current vanishes. In [12] and [2U] this has been considered in the case of two fluids. 
What was found is that isotropization occurs if at least one of the two fluids has a speed of 
sound which is less or equal | the speed of light. It is shown in particular for the case of two 
pressure free fluids in [20] , which can be seen as a singular solution of the Einstein- Vlasov 
system. 

Of course there are many other ways of generalizing results which have been obtained for 
(single) perfect fluids. See for instance [S] and references therein for the inclusion of a 
Maxwell field in the Bianchi I case. In presence of a cosmological constant the results 
of [9] have been generalized even to the Einstein- Vlasov-Maxwell case [11]. A natural 
generalization of the Vlasov equation is the case where the collision term is not zero, i.e. 
the Boltzmann equation. For this case dust-like asymptotics have already been obtained 
in [2T] for an isotropic spacetime with a cosmological constant and [TU] provides a basis for 
a possible extension to the asymptotics in the case of Bianchi I with LRS symmetry. 
Finally we would like to mention that non-diagonal Bianchi I spacetimes are not only 
of interest in the context of cosmology, see for instance a recent work on the so called 
ultra-local limit [3]. 
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